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UNSTEADY HYBRID VORTEX TECHNIQUE FOR TRANSONIC VORTEX FLOWS 

AND FLUTTER APPLICATIONS 

By 

Osama A. Kandil* 

SUMMARY 

This report covers the progress of research work performed under this grant 
in the period of January 1, 1987 to July 31, 1987. During this period of 
seven months, the following papers which report on the accomplishments and 
results have been presented and published (copies of some of these papers 
are attached to this report): 

I. Integral Equation Solution of the Full-Potential Equation for Steady and 

Unsteady Transonic Flows: 

A. "Integral Equation Solution For Transonic And Subsonic Aero- 
dynamics." Kandil, 0. A. and Hu, H., presented at the Third GAMM- 
Seminar on Panel Methods in Mechanics, Kiel, West Germany, January 
16-18, 1987. Published in "Notes on Numerical Fluid Mechanics". B. 

B. "Full-Potential Integral Solution for Transonic Flows With And 
Without Embedded Euler Domains," Kandil, O.A. and Hu, H., AIAA 87- 
1461, Honolulu, Hawaii, June 8-10, 1987. Submitted for publication 
to AIAA Journal. 

C. "Transonic Airfoil Computation Using the Integral Equation With and 

Without Embedded Euler Domains," Kandil, O.A. and Hu, H., 9th 

International Conference on "Boundary Element Methods in 

Engineering, University of Stuttgart, Stuttgart, West Germany, 

♦Professor, Department of Mechanical Engineering and Mechanics, Old Dominion 
University, Norfolk, Virginia 23508 


E. "Influence of Numerical Dissipation on Computational Euler Equations 
for Vortex-Dominated Flows," Kandil, O.A. and Chuang, A.H., AIAA 
Journal, Vol. 25, No. 9, September 1987. 

III. The following abstracts have been submitted for presentations; 

A. "Unsteady Vortex-Dominated Flow Around Maneuvering Wings Over a Wide 
Range of Mach Numbers," Kandil, O.A. and Chuang, A.H., submitted for 
presentation at AIAA 26 ASM, Reno, Nevada, January 11-14, 1988. 

B. "Unsteady Transonic Airfoil Computation Using the Integral Equation 
Solution," Kandil, O.A. and Hu, H., AIAA 27 SDM Conference, 
Williamsburg, Virginia, April 19-21, 1988. 

IV. The following talks will be given at the ICASE, Nasa Langley: 

A. "Ful 1-Potential Integral Solutions for Steady and Unsteady Transonic 
Airfoil computations With and Without Embedded Euler Domains," Hu, 

H. and Kandil, O.A. , August 18, 1987. 

B. "Steady and Unsteady Euler Computations of Various Flow Regimes Past 
Delta Wings,” Chuang, A.H. and Kandil, O.A., August 27, 1987. 

V. The Principal Investigator has chaired and organized the following 

sessions and conferences of the AIAA: 

General Chairman and organizer. Fluid Dynamics Sessions (29 Sessions), 

AIAA 25th Aerospace Sciences Meeting, Reno, Nevada, January 12-15, 

1987. 
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B. Chairman and organizer. Incompressible Navier Stokes II , AIM 8th Compu 
tational Fluid Dynamics Conference, Honolulu, Hawaii, June 9-11, 1987. 

C. Chairman and organizer, Leading-Edge Vortex Flows , AIM 9th Fluid 
Dynamics, Plasma Dynamics and Lasers, Honolulu, Hawaii, June 8-11, 

1987. 

D. Chairman and organizer. Vortex Flow II , AIM 25th Aerospace Sciences 
Meeting, Reno, Nevada, January 12-15, 1987. 
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INTRODUCTION 

Computation of the transonic flow around airfoils and wings 
using finite-difference and finite-volume methods requires 
fine grids and large computational domains around the source 
of disturbance. The outer boundary of the computational 
domain is usually placed at several chord lengths away from 
the inner boundary. Moreover, special treatment is required 
at the outer boundary to approximately satisfy the farfield 
boundary conditions. Using different levels of 

approximations, lnvlscid computational schemes have been 

developed based on the Transonic Small Perturbation (TSP) 
equation (Murman and Cole* ; Edwards, Bland and Siedel^), Full 
Potential (FP) equation (Steger and Lomax' 3 ; Garabedlan and 

Korn 4 ; Jameson^) and Euler equation (Jameson**). These schemes 
require large capacity of computer memory to handle the large 
number of grid points and the associated flow variables. They 
also require large CPU time to obtain converged solutions due 

to the large number of iteration cycles usually it is of 

order a thousand. 

The potential equations can be used for flows with weak 
schocks since the entropy increase and vorticity production 
across these shocks are small. For strong shocks, the 
irrotational and isen tropic flow assumptions are invalid and 
one cannot use the potential flow equations across or behind 
the shocks. For these flows, one has to correct the potential 
equations in order to include the entropy jump across the 
shock wave. Methods of this type have been developed by Hafez 
and Lovell^, Fuglsang and Williams® and Whitlow, Hafez and 
Osher*. Alternatively, one has to use the computationally 


more expensive Euler equations. 

The Integral Equation Formulation (IEF) using the TSP 
(Piers and Sloof 10 ; Tseng and Morino 11 ) or the FP (Kandil and 
Yates^; Oskam*^; Erickson and Strande 1 ^; Sinclair*^; Kandil 
and Hu***) equations represents an alternative to the finite 
difference and finite volume methods to treat transonic flows 
with weak shocks. With the IEF, the farfield boundary condi- 
tions are automatically satisfied and only a small computa- 
tional domain is needed around the source of disturbance. 
Moreover, the accuracy of the method depends on the evaluation 
of integrals rather than derivatives and hence coarse grids 
can be used within the small computational domain. Because of 
these obvious advantages of the methods which are based on the 
IEF, it is highly desirable to fully develop these methods and 
extend them to treat transonic flows over a wide range of Mach 
numbers. 

In this paper, we present a transonic integral equation 
method which is based on the full potential equation, and 
couple the method with embedded Euler domains to treat strong 
shocks. The method is extensively applied to different 
airfoil sections and Mach numbers. The results are compared 
with experimental data and other results which were obtained 
by using finite-difference and finite-volume methods with TSP, 
FP and Euler equations. 

FORMULATION 


Full Potential Equation 

For transonic flows with weak shocks, the dimensionless 
governing equations of the two-dimensional, steady, potential 
flow are given by 


9 + $ 

xx 


yy 


(i) 


G - T (p * ®* + p y V <2> 

P - [1 +^M 2 (1 -®y)] l/rl . (3) 


where <5 is the total velocity potential, p the density, 
the frees tream Mach number, y the ratio of specific heats and 
the subscripts x and y refer to the derivatives. It should be 
noted that G represents the total compressibility in the flow; 
The characteristic parameters are the frees tream velocity 
(U ) and density (p ) and the airfoil chord length (A). 
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The boundary conditions on the airfoil g, away from g and at 
the trailing edge TE are given by 


§ 

• 

H 

O 

on g(x, y) » 0 

(4) 

V$ e 

CD 

away from g 

(5) 

A Cl - 0 


(6) 

P 'TE 



Vg/|vg|, e,, Is 

a unit vector parallel to U 

and 


ACp is the pressure jump. The pressure coefficent is given by 

2 r Y— 1 2 2 2 

c •— h {[1 +V-N. (1 - it - «,)] -1) (7) 

P y M 7 

1 CD 

The formal solution of Eq. (1) in terms of the velocity field 
V® with explicit contribution of the shock surface S is 
obtained as 


v®<*,y) - 1. + ^ i V »> ~ I 7 ' 11 !; 3 ds 

g 6 (x-C) + (y— n) 



where q g and y are the airfoil surface distributions of 
sources and vorticity; respectively, and q^ is the source 
strength of the shock surface. In the shock-capturing scheme, 
the last integral term In Eq. (8) is dropped, since this term 
is Included in the third integral term of the equation. In 
the shock capturing- shock fitting scheme, the last integral 
term corresponding to the shock .surface is retained when shock 
fitting is used. For shock fitting, the following equations 
are used to determine the shock strength q^, properties behind 
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the shock p^, V 2tx* V 2t an< * t * ie or * entat * on the segments 
forming the shock surface 0: 


2 V 


«1. - V 2«> 


In 


Y+l 


(1 - 


M 


In 


2n 


(Y-D Mj n + 2 

(T+l) «L 


In 


M ln > 1 


(9) 


( 10 ) 


V a V 
v 2t It 


( 11 ) 


( T*> M in 

(Y-D M^ n + 2 


( 12 ) 


p » sin [ 


-1 rl.2 sing slnQ . 1 


cos (p-9) 


-41 


1/2 


m: 


(13) 


M, 


r± ni 

MA 2 • "u ■ 7 ®1 • V P 1 2 


(14) 


where the subscripts 1 and 2 refer to the conditions ahead and 
behind the shock, respectively, while n and t refer to the 
normal and tangential directions to the shock surface S, 
respectively, and 9 is the relative direction of the flow 
behind the shock to that ahead of the shock. 


Euler Equations 

For strong shocks, an embedded computational domain is 
constructed around the shock which has been preliminary found 
by the Integral solution with shock capturing only. With the 
boundary and initial conditions found from the Integral 
solution, the unsteady conservative form of the Euler 
equations are solved in this limited domain with psuedo time 
marching. The dimensionless conservative form of these 
equations is given by 


M + M + 2L m 

dt 5x 5y 


0 


(15) 


where the flow vector field q and the flux components E and 
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F are given by 


q 

E 

F 


[p, pu, pv, pe] t 


(16) 

[pu, pu 2 + p, puv , 

puh] t 

(17) 

[pv, puv, pv 2 + p. 

pvh] C 

(18) 


The total energy and enthalpy per unit mass are given by 

6 “ (y+i) p + + y2 ^ 2 * h " e + p/p (19) 


Since we are interested in the steady flow solution only, the 
energy equation [last elements in Eqs. (16)-(18)], which is a 
differential equation, is replaced by the algebraic steady 
form which states that the total enthalpy is constant. Hence, 
the energy equation is replaced by 

p - (p/y) [1/M* + (y-1) (1 - u 2 - v 2 )/2] (20) 

METHOD OF SOLUTION 

Shock-Capturing Shock-Fitting (SCSF) Scheme 

The basic difference between the incompressible Integral 
Equation Solution and the transonic Integral Equation Solution 
are the additional third and fourth integral terras of Eq. (8). 
In the shock capturing part of the scheme, the fourth integral 
terra is dropped while in the shock fitting part of the scheme 
this term is taken into account. 

The SCSF-scheme is an iterative scheme due to the 
nonlinearity of the third and fourth integral terms. The 
iterative scheme is described below: 

Neglecting the fourth integral term and setting G ■ 0, a 
standard panel computation is used to obtain q g and/or y . 
The q and y are piecewise linear distributions on each 
surface panel ^and they are defined in terms of their nodal 
values. Initial values of G are calculated at the centroids 
of the field elements by using the linear compressibility 

relation G - M 2 u . where u„ Is the x-derivative of the x- 

CD X X 

component of the velocity. The centroidal value G represents 
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the G value for the field element. Equation (8) Is then used 
to enforce the surface boundary conditions, Eqs. (5) and (7) 
to find new q g and/or y g . The density p and nonlinear 
compressibility G are calculated by using Eqs. (3) and (2). A 
type finite-difference expression is used to calculate p and 
p depending on the type of the centroidal point-subsonic or 
supersonic. Once the G values are obtained, the surface 
boundary conditions are satisfied again. The iterative 
procedure is continued until the shock location is fixed. 
This is the shock capturing part of the scheme. 

Shock panels are then introduced at the shock location, 
the fourth itegral term of Eq. (8) is now taken into account, 
Eqs. (9) and (13) are used to calculate q 3 and p and Eqs. 

(10)-(12) are used to cross the shcok panel. The iterative 
procedure is continued as before with the exception of dealing 
with the shock panels as explained. Convergence is achieved 
once the surface pressure converges. This is the shock 
fitting part of the scheme. 

Integral Equation With Embedded Euler (IEEE) Scheme 
In this scheme, the shock capturing part of the SCSF-scheme is 
used to locate the shock. Once the shock is captured, a fine 
grid is constructed within a small computational region around 
the shock where a finite-volume Euler scheme is used. The 
basic f ini te-volume equation is obtained by integrating Eq. 
(15) over x and y to obtain 


// ^ dA + <J) (E dy + F dx) - 0 (21) 

Equation (21) is then applied to each cell of the embedded 
grid of the Euler domain. The resulting difference equation 
is 

4 

(§J) A A + Z (E A y + F A x ) ■ 0 (22) 

5t u 13 t-i r r 

where AA is the cell area, r refers to the cell-side number 
and the integer subscripts i, j refer to the centroidal 
values. The Euler solver is a central-difference finite- 
volume method which uses four-stage Runge-Kutta time stepping 
with explicit second- and fourth-order dissipation terms. The 
details of this solver are given in reference [17]. 

The boundary and initial conditions for the Euler domain 
are obtained from the Integral equation solution which is 
interpolated on the Euler domain grid. The Euler solver is 
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then used to capture the shock and calculate the flow vector 
field q. It should be emphasized here that the downstream 
boundary condition must be_ updated while Euler calculations 
are excuted. Fixing the q values of the Euler domain, the 
IE calculation is used to update the boundary conditions for 
the Euler domain. The iterative procedure is repeated until 
convergence is achieved. 

COMPUTATIONAL EXAMPLES 

In this section, we present applications of the SCSF- and 
IEEE-scheraes to the NACA 0012 and 64A010A. According to the 
convergence study using different sizes of the IE 
computational domain, which was presented by the authors 

(Kandll and Hu 1 ^), a computational domain of 2 x 1.5, in the x 
and y directions, has been used around the airfoil in all the 
following applications. A rectangular grid of 64 x 60 has 
been used for the IE computation. The grid is clustered in 
the leading edge, the shock region and near the airfoil 
surface. Since the third Integral terra of Eq. (8) is 
computationally expensive, its computation with constant G 
distribution has been restricted to the nearfield computation. 
For the farfield computation, this term is replaced by the 
equivalent lumped source term at its centroid. With 
sufficient accuracy, it has been computationally determined 
that the nearfield distance from. the centroid is < 0.5. 

Figure 1 shows the results of the SCSF-scheme for NACA 
0012, M * 0.8 and a ■ 0° , along with comparisons with the 

IQ 

computational results of Garabedian, Korn and Jameson , and 
the experimental data taken from reference 19. The SCSF 
scheme took 12 iteration cycles of shock capturing (SC) and 13 
cycles of shock fitting (SF) to achieve convergence. 

Figure 2 shows the results of the IEEE-scheme for the 
same case along with a comparison with the computational 
results of Jameson** who also used the finite-volume Euler 
scheme with four-stage Runge Kutta time stepping. In the 
present IEEE-scheme, the embedded Euler domain has a size of 
0.5 x Q.6 around the shock region with a grid of 25 x 30. 
This case took 10 iteration cycles of SC, 250 time cycles of 
Euler iterations to achieve a residual error of 10~^ and 5 IE 
cycles to update the Euler domain boundary conditions. 

Figures 3 and 4 show the results of the SCSF- and IEEE- 
schemes for NACA 64A010A, M o a 0.796, a ** 0° along with 
comparisons with the computational results of Edwards, Bland 
and Seidel 2 who used the TSP equation, and the experimental 
data taken from reference 2.® With the SCSF-scheme, the 
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numbers of SC and SF iteration cycles to achieve convergence 
are the same as those of the case presented in Fig. 1. With 
the IEEE-scheme, the embedded Euler domain has a size of 0.7 x 
0.6 with a grid size of 35 x 30. This case. Fig. 4, took 10 
iteration cycles of SC, 130 time cycles of Euler iterations to 
achieve a residual error of 10"^ and 3 IE cycles to update the 
Euler domain boundary conditions. 

Figures 5 and 6 show the results of the SCSF- and IEEE- 
schemes for the lifting case of NACA 0012, « 0.75 and 

a » 2° along with the computational results of Steger and 
Loraax^ , and the experimental data taken from the same 
reference. The size of the grids and the number of iteration 
cycles used to achieve convergence are the same as those of 
the cases given in Fig. 1 and 2. 

Figure 7 shows the results of the IEEE for NACA 0012, 
■ 0.812 and cc ■ 0° along with the experimental data of 
reference 18. In Fig. 8, the results of the IEEE for NACA 
0012, ■ 0.82 and a ■ 0° are shown along with the three- 

dimensional solution at the wing root chord of Tseng and 

Morino^, who use the IE for the TSP, and the experimental 
results of reference 20. The size of the embedded Euler 

domain for these cases is 0.8 x 0.8 end the computational grid 
is 40 x 40. 

CONCLUDING REMARKS 

Two transonic computational schemes which are based on the 
Integral Equation Formulation of the full potential equation 
have been presented. The first scheme is a Shock Capturing- 
Shock Fitting (SCSF) scheme which uses the full potential 

equation throughout with the exception of the shock wave where 
the Rankine-Hugoniot relations are used to cross and fit the 
shock. The second scheme is .an Integral Equation with 

Embedded Euler (IEEE) scheme which uses the full potential 

equation with an embedded region where Euler equations are 
used. The two schemes are applied to several transonic 

airfoil flows and the results have been compared with numerous 
computational results and experimental data. The two schemes 
are nevertheless efficient as compared to the other existing 
schemes which use finite-difference or finite-volume methods 
throughout large computational domains with fine grids. The 
SCSF-scheme is restricted to flows with weak shocks, while the 
IEEE-scheme can handle strong shocks. Currently, the IEEE 
scheme is applied to other transonic flows with strong shocks 
as well as to unsteady pitching oscillations. 
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OUTLINE OF THE TALK 


1- Motivation and Objectives 
2* Formulation 

- Classical and Zero-Total-Pressure-Loss Sets 

- Supersonic Conical Flow Equations 

- Relative Motion in a Rotating Frame of Reference 


3. 

4. 


Highlights of Method of Solution 


Applications : 

- Conical Flow, Sharp-Edged Wings (Classical and ZT PL Sets) 

- Conical Flow, Round-Edged Wings (Classical and ZTPL Sets) 

- Three-Dimensional Flows; Transonic and Low-Speed Flows 

- Uniform Rolling in a Conical Flow 

- Rolling Oscillation in a Locally Conical Flow 


Ex* — p ^ 5 (0 


fs/u*»«r;c mS 


5« Concluding Remarks 



1 Classical 
Vortex 

2 Separation 
Bubble with 
No Shock 

3 No Shock/ 

No Separation 

4 Shock with no 
Separation 

5 Shock-Induced 
Separation 

6 Separation 
Bubble with 
Shock 

7 Vortex with 
Shock 


Figure 1 . Miller and Wood 1 Classification Diagram. 



CLASSICAL EULER EQUATIONS 


• Conservation Form of Euler Equations in a Space-Fixed Frame of Reference 


aq * at * a? . aG n 
U + + Ty + m 0 


(i) 


5 ■ Cp. pu, pv, pw, pe] T 


( 2 ) 


t » [pu, pu 2 ♦ p, puv, puw, puh]* 


( 3 ) 


f • [pv, puv, pv 2 ♦ p, pvw, pvhj 1 


( 4 ) 


5 • [pw, puw, pvw, pw 2 ♦ p, pwh ] 1 


( 5 ) 


e • p/ptr-D ♦ ( u 2 ♦ v 2 ♦ w 2 )/2 


( 6 ) 


h * e + p/p 


( 7 ) 


Conical Variarles 


SUPERSONIC CONICAL FLOW EQUATIONS 


£ • X , n « y/x , C ■ 2 /x 


( 6 ) 


Conical Flow Equations 


si ♦ S * H ♦ 2 e * o 


( 7 ) 


WHERE 


P« F * T) [ 


( 8 ) 


g - G - C E 


( 9 ) 


Zero-Total-Pressure-Loss Euler Equations 


• Replace the energy equation by either 
one of the isentropic gas equations 

p/p Y * const, or + V *(psv) = 0 


Set (1): • Replace the x-momentum equation (second 

elements in the vectors given in Eqs. 
(2) ‘and (3)) by the steady energy 
equation (total constant enthalpy) 

h * const * , r- + \ (u 2 + v 2 + w 2 ) 
TpTTp 2 


Set (2): • Replace the continuity equation (first 

elements in the vectors given in Eqs. 
(2) and (3)) by the steady energy 
equation given in Set (1) 


EXPLANATION OF TOTAL-PRESSURE CHANGE 
FOR CLASSICAL AND ZTPL SETS OF EULER EQUATIONS 

Differential Euler Equations 
Crocco's Theorem 


I ? s * * * ? ♦ ♦ n 

Definition of Entropy Change 


A S * R In - — + C In 

P T P o 


(a) Classicai Set 


Steady flow ^ * o, h - const 


and 


TVS-wxY. AS-Rln 


( 1 ) 


< 2 ) 


For a free-sheet Z is parallel to V! 7 S 3 0 - Pj * Py ♦ Zero-Total-Pressure Loss 
(p) Zero-total-Pressure-Loss Set (Shock-Free and Weak Shocks) 

h * const, 7 S ■ 0 

and must pe parallel to V, Pj * Pj ♦ Zero-Total-Pressure Loss 


Cqhputatqnal Euler Equations 
Crocco's Theorem 


T a s ■ « x v + ^ ♦ vh . I v • ? Viscous-Form of the Equation 


(l) 


Definition of Entropy Change 


a $ » r in * — + cx 


T P n V 


(2) 


(a) Classical Set 

For STEADY FLOW - 0, h » const AND 


TVS»uxV + Numerical Dissipation, A S * R Xn 


Even if w is parallel to V, v S * 0 ♦ Pt * P T ♦ Non-Zero TPL 

* m 1 

(B) Zero-Total- Pressure-Loss Set (Shock-Free and Weak Shocks) 

' h ■ const, A S * 0 

0 » u x V + Numerical Dissipation, P^ ■ P^ * Zero TPL 

i 


CLASSICAL EULER EQUATIONS FOR THE RELATIVE MOTION IN A ROTATING 
FRAME OF REFERENCE 

• The conservation form of the classical Euler eouations for the absolute motion 

OF THE FLOW IN A SPACE-FIXED FRAME OF REFERENCE 


♦ v • (p v) • o (i) 

? ( P ♦ 7 • (p V V + p T) - 0 (2) 

o t 

H g -l l + V • (p h V) « 0 . (3) 

V 2 

e ■ p/p ( y- 1 ) ♦ y- 

h ■ e ♦ p/p ( 5 ) 


I 



To express these equations in terms of a 
moving frame of reference, we use the following 
relations of the substantial and .local 
derivatives for a scalar "a" and a vector "A": 


Da . D'a 
Ft ' FT" 


aa a d'a 
at at r 



7a 


DA 

Ft 



wxA 


3A a'A 
at “ FT 


V t • 7A + uxA 


(6. a) 


(6.b) 


(7. a) 


(7.b) 


• The transformation velocity V is a function of 

this moving frame of reference translation and 
rotation 


V = V - V = V + wxr (8) 

t r o 


Restricting the motion of the frame of reference to the 
rotational motion, 

OV 

V o = 0 and = q, 

the equations of relative motion in the rotating frame of reference 

+ 7 • (p V ) » 0 
Ot r 

a'(pv) 

§p + 7 * [p V f V p + p I] , _ pfuxf + 2toxV + 


S' (p e ) 

— ap— + 7 • [p h. v J 


where 


v r 1 - - 2 

e r * p{ Y - n * T ~ 1 l uxr i = e - V • (ujxr) 


(9) 


r + o«(wxr)] (10) 

r v r J s - pCV. • (ix?) + (wxr) • (*xF)] (H) 


(12) 


2 

h r ■ + T ■ 7 M ■ h - V • (uxF) 


(13) 




The 

the 


abstract conservative form 
rotating coordinates (x' t y'. 


a 'V, s ' J r »' p r «' 5 r 

-5T^ * -XT' * Sr~ * 


of the relative 
z') is given by 

S 


motion 


in 


terms of 


(14) 


where 


q r - Cp, pu r , pv r . pw r , 

— 2 

£ r * Cpu r . P u r ♦ p, P u r v r . pu r w r . pty.h r J 1 

* Cpv r * pu r v r* pv r + p * pv r V pv r h r jt 
g r * Cpw r . pu r w r . P v r w r . pw p 2 ♦ p, pWr h r ] t 

S * [0. 0. p(J z ♦ anv r + J y] , _ p{ - y + a-v . w 2 

r 

• Since only the rolling motion is solved 

term S has been written for Z = « e , , and u = i e . 

• X ^ 




p(-V r W2 + W f 

the source 


(15) 

(16) 

(17) 

(18) 

• *2 • p - 1 

ojy + ujuy + (JW2 ) J 

(19) 


H16HUGHTS OF METHOD OF SOLUTION 


1. We use the Central-Difference Finite-Volume Scheme with Four-Stage Runge Kutta 

TIME STEPPING AND EXPLICIT SECOND* AND FOURTH-ORDER DISSIPATION TERMS- 

2- For steady flows, local-time stepping is used, and for unsteady flows minimum 

GLOBAL TIME STEPPING IS USED- 

3* A THREE-DIMENSIONAL COMPUTER PROGRAM IS USED TO SOLVE FOR: 

- Conical flows (using 3 conical planes, we enforce the absolute flow vector to be 

EQUAL ON THESE PLANES) 

- Direct solution of the three-dimensional flow problem- 

u. Depending on the problem under consideration, different initial conditions are 
USED- 

5- Depending on the problem ’under consideration, different surface, farfied and 

SYMMETRY COND I T I ONS ARE USED- FOR SUPERSONIC FLOWS, THE OUTER ROW SHOCK IS 
CAPTURED AS PART OF THE SOLUTION- 


NUMERICAL EXAMPLES (1) 


SHARP-EDGED H1NGS (CLASSICAL EULER EOS- ( ZERO-TOTAL-PRESSURE-LOSS SETS) 

ROUND-EDGED MINGS (CLASSICAL EURER EOS- * ZERO-TOTAL-PRESSURE-LOSS SETS) 

- NUMERICAL BOUNDARY CONDITION (COARSE AND FINE 6RIDS) 

- CLOSED FORM BOUNDARY CONDITION (COARSE AND FINE GRIDS) 


THREE-DIMENSIONAL TRANSONIC AND SUBSONIC FLOWS 


original, page ® 

OF POOR QUALITY 


Figure i . Standard Euler Set, Sharp-edged Wing, 64X64 Cell, 

^<»* 2.0, a»10®, (1*70°, £ 2 * 0 . 12 , £^*0.005, 1. Surface pressure, 2 
Crossflow Mach number, 3. Crossflow Velocity. 








Figure 2. Zero-Total -Pressure-Loss Euler, Set (1), Sharp-edged Wing, 64x64 
Cell, M - *2.0, a*10°, p*70°, c 2 *0.12, e 4 *0. 005,1. Surface Pressure 

2. Crossflow Mach number, 3. Crossflow Velocity. 

i in 1 rft 1 KM/ 2 / I Mil U m 
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Invite Id Shock. 

Induced Separation 




Details ef Laadtng.cdga Flew (r, f . 2.3) 



